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ABSTRACT 


Composite materials have large number of uncertain parameters associated with its 
production. It is not possible to control these completely. This results in variations in 
the properties of composite materials. For accurate modeling, the material properties are 
considered as random \-ariables. The influence of randomness in material properties on the 
deflection behaviour of T300/52208 graphite/epoxy composite laminates has been analysed 
in the present stud}' with the help of Monte Carlo Simulation and FEM. The input material 
properties in a lot of samples are assumed to vary according to normal distribution and 
the effect of this is investigated on deflection of plates with different boundary conditions, 
aspect ratios, end conditions and ply-orientation. The study has also been extended to 
plates with cut-outs. It has been found that the normalised displacement covariance of the 
plate reacts to the input randomness in the same way for various aspect ratios, thickness 
ratios, ply-orientations and support condition. This results in the development of a design 
curve with which the deflection statistics of a different plates can be predicted with a fair 
degree of accuracy. 
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Chapter 1 


Introduction 


Composite materials are known to have many advantages over conventional materials. 
Further these allow the designer the flexibility of tailoring them to meet specific needs. 
This makes composites a contender for the material of the future, which is also evident 
from the fact that more and more composite materials are replacing conventional materials 
in engineering and industrial applications. The mechanical and physical properties of the 
raw materials fiber, resin etc. as provided by the manufacturers are generally average 
properties subject to specific production conditions. Because of inherent variations in 
these conditions from time to time a dispersion in the properties of raw materials is always 
present. 

Fiber and resin are combined in a suitable proportions to form a lamina. Factors such 
as temperature, pressure, moisture, fiber spacing, voids, impurities curing time etc. enter 
along with the variations in raw material properties. These manifest themselves as the 
variation in lamina properties. The lamina are stacked together to meet the requirements of 
the designer to form laminates, the variations in the lamina properties results in variations 
in the reduced stiffness coefficients of the laminates. 

As shown in Fig. 1.1 the uncertainty in the material properties gets introduced at various 
levels through different parameters associated with the composite. Further it is well under- 
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Table 1.1: Variation of reported stiffness properties for a widely used Graphite/Epoxy 
system , from five sources (Sources: Major Airframe Company Reports). 


Elastic Constants(10®psi) 

1 

2* 

3* 

4 

5 

Max. Diff.% 

Long, tensile modulus 

20.8 

18.1 

21 

20.6 

18.5 

16 

Long, compr. modulus 

18.6 

14.5 

21 

19.8 

18.5 

45 

Trans, tensile, modulus 

1.9 

1.8 

1.7 

1.3 

1.6 

46 

In-plane shear modulus 

0.85 

- 

0.65 

0.8 

0.65 

31 

Poisson’s ratio 

0.30 

- 


0.32 

0.25 ' 

28 


-Different divisions of the same company 


stood that it is not possible to control all the parameters involved in the manufacturing of 
composites. Tables 1.1 and 1.2 [1] illustrate the variation in strength and stiffness proper- 
ties of composite laminates, observed in a molding lot. The variations in the strength and 
stiffness properties introduces a factor of uncertainty in the behavior of structures made 
up of the composite materials. Accurate predictions of the system behavior of composites 
in the presence of the uncertainties favors a probabilistic analysis approach for composites 
by modeling their properties aa random variables (RVs). 

The uncertainty in the material strength and stiffness of composites has caught the 
attention of researchers. Ibrahim [2] reviewed a number of topics pertaining to struc- 
tural dynamics with parameter uncertainties. Nakagiri et al. [3] have studied simply sup- 
ported(SS) graphite/epoxy plates with Stochastic Finite Element Method (SFEM) taking 
fiber orientation, layer thickness and layer number as random variables, and found that 
the overall stiffness of FRP laminated plates is found out to be largely dependant on fiber 
orientation. Leissa and Martin [4] have analyzed the vibration and buckling of rectangular 
composite plates, and have established that variation in fiber spacing or redistribution of 
fibers tend to increase the buckling load by 38 percent and the fundamental frequency by 
21 percent. Englestad and Reddy [5] have studied metal matrix composites (MMC) based 
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Table 1.2; Variation of reported stiffness properties for a widely used Graphite/Epoxy 
system, from five sources (Sources: Major Airframe Company Reports). 


Strength Properties(10®psf) 

1 

2* 

B 


5 

Max. Diff. % 

Long, tension 

274 

190 




67 

Long, compression 

280 

126 

180 

126 

162 

122 

Trans, tension 

9.5 

5.2 

8.0 

5.4 

6.0 

83 

Trans compression 

39 

- 

30 

21 

25 

86 

In-plane shear 

17.3 

- 

12 


- 

106 

Interlaminar shear 

- 






90 


^-Different divisions of the same company 


on probabilistic micromechanics nonlinear analysis procedure of MMC. They have used 
Monte Carlo Simulation (MCS), and different probabilistic distributions to incorporate 
the uncertainty of constituent level properties. Roy and de Wilde in first part of their 
work [6] have established a procedure to obtain the best fit for each material property. 
Further they have studied the behavior of perforated plate and determined the probability 
distribution of the response. The input variables have been simulated using MCS and 
SFEM has been used for the analysis. Salim et. al [7, 8, 9] studied the statistical response 
of the system considering the material properties such as longitudinal modulus, Poisson’s 
ratio, transverse elastic modulus, etc., as independant random variables on [90°/0°/90°] 
and [0°], rectangular, SS and clamped(CC) plates subject to static loads. Salim et 0^9] 
have shown that there is change in the rate and the order of the natural frequency with 
SD of input RV and change in boundary condition of the plate. 

One can see that the uncertainty in composite material properties enters at the raw ma- 
terials level [4, 5|. In references [1, 3, 6] and [7] the dispersion in the behavior of composite 
structure is assumed to be due to the variation in lamina properties. These studies also 
takes care of the variation in the properties of the raw materials. Salim et a^7] used the 













Manufacturing Tolerance 



Figure 1.1: Development of uncertainty in the material property of laminated composites. 
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perturbation technique with Rayleigh-Ritz formulation to analyse the bending, buckling 
and vibration of composite plates. With the Rayleigh-Ritz method, different formulations 
are needed to handle different problems. That is the formulations have to be changed for 
different end conditions, different loadings, etc. This problem can be over come by using 
FEM for deflection analysis. The perturbation technique though gives acceptable results, 
the method fails for larger dispersions in the input random variables . The SFEM method 
can handle larger uncertainties in the input parameters. However, it is computationally 
expensive and being a specialised technique its other applications are limited. This has 
motivated the use of a conventional FEM in the present study. Though this method is 
also computationally intensive, the analysis can be done with commercially available FEM 
modules. Most of the earlier studies have been performed where material properties vary 
independantly. In a practical problem the material properties vary randomly and simulta- 
neously. The present study attempts to cover this aspect also. 

The variation in the properties are present in a given lot of plates as well as there are 
variation of the individual properties over the domain of the plate. The present study 
focusses on the variation in the lot, variation in property in a given lot which is assumed 
to conform with the normal distribution. In the present study the stiffness properties of 
lamina such as Et, uu, Gu-, Gtz and G^ are considered to be the basic random variables. 
Deflection is one of the important behavioral parameters of a structure and this may be 
a constraint on the design for certain applications. In the present analysis the bending of 
composite laminates due to external loading has been studied with the above mentioned 
material properties modelled as random variables. Effect of the dispersion in the input 
material properties is studied on the flexural response of the plate. MCS with FEM is 
adopted to obtain the deflection response of the laminated plates and its statistics with 
different input characteristics and cut out shapes. 


Chapter 2 


Formulation for static analysis 


In this section, the details of finite element analysis for plate bending is being presented. 
A suitable displacement model is adopted and from these the strain displacement relations 
are developed. Since this study is concerned with composite plates, the stress-strain rela- 
tionship for an orthotropic lamina is presented. These are utilized to obtain the laminate 
constituent equation. Finally the governing equation for the flexural response of laminated 
anisotropic plates is presented. 

2.1 Displacement field 

The rectangular laminated plate analyzed is represented in Fig.2.1. u, v and w are the 
displacements along x, y and z axes respectively, uq, vq and Wq represent the mid-plane 
displacements, and ipy are the rotations. Thus the nodal displacements model has five 
degree of freedom namely, «o, Uq, wq, ipx, iiy 

A displacement field of the following form is used for the present study. 

u(x, y, z) - «q{x, y) + fi{z)'tpxix, y) + h{z)dwQldx 

v{x,y,z) = VQ{x,y) + fi{z)ipy{x,y) + f2{z)dwo/dy (2.1) 

wix, y, z) = wq{x, y) 
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where, 

/l(z) = CiZ-C2Z^ 

/a ( 2 ) = -CiZ^ (2.2) 


The displacement field can be transformed either to first order shear deformation theory 
(FSDT) model [ 10 ] or higher order shear deformation theory {HSDT)model by changing 
the values of the constants Ci, C 2 and C 4 [ 11 , 12 ]. 
with, 


Cl = 1 ;C 2 = €4 = 0 for FSDT model 
Cl = 1 ;C 2 = C 4 = -:^ for HSDT model 

In the displacement model given above, u and v are dependent on the derivative of the out 
of plane displacement w. The strain vector would have a second order derivative of the out 
of plane displacement hence for the FEM analysis its necessary to employ an element with 
a C^ continuity . To overcome the complexity and difficulty involved with making a choice 
of C^ continuity element, the derivative of the out of plane displacement are themselves 
considered as separate DOFs 6x, dy. Thus the 5 DOF with C^ continuity is transformed 
into a 7 DOF C° continuity. 

Thus the displacement model adopted in the present study is given by the following : 



u{x, y, z) = Uo(x, y) + fi{z)‘ipx{x, y) + f2{z)9x{x, y) 

v{x,y,z) = vo{x,y) + fi{z)i)y{x,y)-{- f2{z)9y{x,y) (2.4) 

w{x,y,z) = wo{x,y) 


where, 9^ = dwoldx and 9y = dwo/dy. 
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2.2 Strain-Displacement relations 

The strain displacement relations are obtained using small deformation theory [13]. 




{ej = 


Ixy 


^yz 

Ifxz 


— {^2} + 


where, ej is the bending strain vector and e, being the shear strain vector. 
Where, 


{^1} = 


{^2} 


eii 
€21 ^ = i 


{«i} = { 




€41 

^51 

/ \ 
Kli 

«21 

K61 


Uq^x 

^0,y "h ^0,1 

% 

'<Px 


Ci{ 


+ < 


u;o,y 

Wo, I 


= Ci 


V / 


'^x,x 

'^y,y 

'^x,y + '4^y,x 



«;i2 


'^x,x 


/ x 

{/C2} = < 

K22 

>=-C2< 


► - <74 < 



'«62 

V . J 


+ ’0y,x^ 


^x,t/ + 

V . 


{«3} = < 

[«4ll 
\ 1 

^ = - 3 C 2 ^ 

\^y] 

1 j 



[«5lJ 


[^xj 



(2.5) 


( 2 . 6 ) 
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2.3 Stress- Strain relations for an orthotropic lamina 

The constituent equation for k*'^ lamina of the laminate can be written as [13] 



k 

r 




•n 

k 

'' \ 

cri 


Qu 

Qi 2 

0 

0 

0 


ez 



Qi 2 

Q22 

0 

0 

0 


ez 

' Tit 

> = 

0 

0 

Qee 

0 

0 

< 

lit ^ 

Ttz 


0 

0 

0 

Qu 

0 


ytz 

Tlz 

J 


0 

0 

0 

0 

Q 55 




where ‘1’ and ‘t’ represent longitudinal and transverse directions and [Qij]$ are the stiffness 
coefficients of the lamina. These are related to matrix properties as follows: 


Qn 

= El! (1 - i/itVti) 

Q22 

= EtHl-uitUu) 

Qi 2 

= vitEtj {1 - vitvti) 

Qee 

= Git 

Qu 

= Gtz 

Q55 

= Giz 


The transformed stress-strain relations with respect to laminate co-ordinate system are, 


N 

(fc) 

Qn 

Qi2 

Qi6 

0 

0 

(k) 

f X 

: % 


Qi2 

Q 22 

Q 26 

0 

0 




> = 

Qu 

Q 26 

Qee 

0 

0 

< 

7xy ^ 

'^yz 


0 

0 

0 

Qu 

Q 45 


lyz 



0 

0 

0 

Q 45 

Qs5 


Txz 

V. / 


(2.9) 
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Here [Qij] are the transformed reduced stiffness coefficients for the Kth lamina. 

Qn ~ QnC^ + 2 (Qi 2 + 2Q^q)C^S^ + Q 22 S‘^ 

Q 22 — QnS‘^ + 2 {Qi 2 + 2Qqq)C^S‘^ + Q 22 C‘^ 

Qi 2 — (Qn + Qi2 ~ 4Q66)C^S'^ + Qi2(C^ + 5'^) 

Qi 6 — (Qn ~ Qi 2 — 2Q66}C'^'S^ + (Qi 2 — Q 22 + 2Q66)CS^ 

Q 26 — (Qn ~ Qi 2 — ‘^Q66)QS^ + (Qi 2 — Q 22 + 2Qse)SC^ (2.10) 

Qee “ (^11 “f Q22 ~ 2Qi2 — 2Qs6)C^S^ + Qg6(S^ + C^) 

Q 44 — Q44Q^ + Qss'S'^ 

Q 55 — * 3445 ^ + <355^'^ 

^45 == (*344 — Q55)CS 

where C = cos 6 and S = sin 6. 


2.4 Laminate constituent equations. 


The stress and strains are function of all the three co-ordinates x, y and z respetively. To 
reduce the complexity of the number of dimensions, statically equivalent force and moment 
system or the concept of stress resultants is used. These are obtained by integrating the 
respective stresses along the laminate thickness direction [13]. 


Nr Mr 


Ny My Py 
Nxy Mxy Pxy 


Qy ^ 
Qx Px 


- = eT ■ 

f' \ 

(^x 

/ 

Txy 

K ✓ 

■ = 2 C 1 

/ 



ik) 


ay\ (l.Z,Z^)dZ 


(k) 

^ (l,Z^)dZ 


( 2 . 11 ) 


Where, NL gives the number of layers and the resultant force and moment equation can 
be written as; {iV} = [I?]{e} 
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where, {A^} = { Ny Nxy My M^y Px Py Pxy Qy Qx Ry Rx} 
and {?} = { eii 621 cei «2i ><•12 1^22 «^62 €41 €51 kh K51} 

where . . . {PxPyPxy} and {RyRx} are higher order moments and shear stress resultants, 
and, 


[AI] [B] [£;] 


[D] = 


[B] [Dl] [FI] 
[E] [FI] [H] 
0 0 0 
0 0 0 


0 0 

0 0 

0 0 

[A2] m 

[D 2 ] [F 2 ] 


The plate stiffness matrices AUj, A 2 ij, Bij, . . . Hij are given by, 


( 2 . 12 ) 


(.41jj Bij Dlij Eij Flij Hij) 


NL .z^ 

X: / Q!^^(l,Z,Z^Z^Z^Z®)dZ,^,j = 1,2,6 

fc=i 


U2ij D2ij F2ij ) 


NL ^Zk 

5 : 4 Qlf(l,Z 2 ,ZVZ,M = 4,5 

Jfc=:l 


Ij 




( 2 . 13 ) 


2.5 Finite element modeling 


The displacement model can be represented concisely by the following: 


{«} = [N]{6} 


where, ()' refers to transpose of the vector/matrix. 

{u} = {u V w}" 


K 


100 0 /2(z) 0 fi{z) 

0 1 0 f2{z) 0 fi{z) 0 

0 0 1 0 0 0 0 


( 2 . 14 ) 


( 2 . 15 ) 


( 2 . 16 ) 
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{5} = {uo Vo Wo 6y 9x ipy TpxV (2-17) 

{5} is the ‘ displacement vector An isoparametric element is employed for the finite 

element modeling. For this type of element the displacement vector and the geometry are 
represented by the same interpolating functions. 

NN 

= (2.1S) 

i=l 

NN 

x = Y,NiXi (2.19) 

i=l 


NN 

y = P-20) 

t=l 

Here NN is the number of nodes per element. Ni is the shape function for the ith node 
and {(5i} is the displacement vector for the ith node. The shape function depends on the 
type of element employed. The shape functions for a nine noded isoparametric Lagrangian 
element is given below 



Ni = 1^(1 -CMi -^) 

n 2 = -\ii-eHi-v) 

Nz = + 

Ns = (2.21) 

N, = + 

Nr = -\m-Oriil + v) 

N& = ^(1 - + 7?) 

Nd = + 0^(1 + n) 
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The strain vectors given in Equation 2.5 and 2.11 may be written as, 


= {Lm 

{ea} = [■^'2]{<5} 

{ ki } = [Ls]{S} ( 2 . 22 ) 

{«2} = 

{fts} = [T5]{(5} 
and{e} = [i?]{(5} 


In the above equations, [Li], . . . [L5] and [£] are matrices of differential operators 


[Li] 

\L2] 


[Lz] 




£ 0 0 0 0 0 0 

0 1- 0 0 0 0 0 

oy 

^ 0 0 0 0 0 

0 0 0 0 Cl 0 

0 0^ 0 0 0 Ci_ 

0 0 0 0 0 0 Ci^ 


0 

0 

0 

0 

0 

Cl 

a. 

dy 

0 



0 

0 

0 

0 

0 

n ^ 

Cl 

a 

dy_ 




0 

0 

0 

0 

C^Tx 


0 

n ® 

— 

0 

0 

0 

^4 

dy 

0 


0 


0 

0 

0 

n ^ 

C4 

dy 

C 

d 

^ dx 



{L 5 ] = 


-3 


0 0 0 ^4 

0 0 0 0 


0 

a 


C 2 

0 


0 

C2 


(2.23) 
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The assembled differential operator is written as: 


A 

dx 

0 

0 

0 

0 

0 

0 

0 

A. 

dy 

0 

0 

0 

0 

0 

j. 

dy 

d_ 

dx 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

0 

c ^ 


0 

0 

0 

0 

-c^i-x 

0 

-C^Tx 

0 

0 

0 


0 


0 

0 

0 

0 




^ d 

0 

0 

j. 

0 

0 

Cl 

0 

0 

0 

d_ 

dx 

0 

0 

0 

< 7 i 

0 

0 

0 

- 3(74 

0 

-3C2 

0 

0 

0 

0 

0 

-3C4 

0 

-3C2 


The strain displacement matrix [B] for the node can be written as, 

[Bi] = [£]Ni (2.25) 

and the element stiffness matrix is, 

{B\[DmdA (2.26) 

Transforming from x, y co-ordinate system to natural co-ordinate system ^, 77 , the 
stiffness matrix is represented as 

Cl Cl 


(2.27) 
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Here [Bi]' is the transpose of matrix [Bj] and, 


iVi.x 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

Ni^y 

Ni,r 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C'iiVi,x 

0 

0 

0 

0 

0 

ClNi^y 

0 

0 

0 

0 

0 

0 

CiNi,, 

ClNl^y 

0 

0 

0 

0 


0 

1 

H 

0 

0 

0 

-C,Ni,y 

0 

-C'aM.y 

0 

0 

0 

0 

-C4Ni^ 

-C^Ni^y 


-C2Ni,y 

0 

0 

Ni,y 

0 

0 

Cl 

0 

0 

0 

Ni,. 

0 

0 

0 

Cl 

0 

0 

0 

-SCiNi 

0 

-ZC2Ni 

0 

0 

0 

0 

0 

-ZC^Ni 

0 

-ZC^Ni 


( 2 . 28 ) 


[J] = 


dx dy 
a? dr} 

dx dy 
-dr) di_ 


9j:ZNrXr 

df} dr} 


When numerical integration is adopted, [k]ij becomes. 


W« “EE W,W,{B.]'lDmdet[J] 

p=l 9=1 



( 2 . 29 ) 


( 2 . 30 ) 


Wp, ir,j are weights used in Gaussian quadrature. 
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2.6 Strain energy of the laminate 

The total strain energy is the sum of strain energy of individual elements 

NE 

U = 

e=l 
NB . 

= EoA (2.31) 

e=l ^ •'•4(e) 

where NE is the number of finite elements into which the problem domain has been divided. 
Substituting the above equation in equation. 2.23, we get 

^ {{^}' 1^1'Plli:] {i}) dA (2.32) 

Substituting for 5 from equation 2.18 in equation 2.32, it may be put as 

= I { (E {•Si}' P) {•5i}) } dA (2.33) 

Further using Eqn.2.25 in Eqn.2.33 yields, 

= 5 ({d}' [Bmm {d}) dA (2.34) 

2.7 Potential energy of the externsd load 

A plate can be subjected to many types of loading such as the in-plane loading, moment 
loading, shear loading and out of plane loading. For the present analysis two types of out 
of plane loading is considered namely the sinusoidally varying and the uniform loading. 
The potential developed due to the external lateral load is given by 

NE 

V = (2.35) 

C=1 

Where, the element load vector is. 




(2.36) 
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Here in this study only the transverse load is applied the loads corresponding to all other 
DOF is kept zero and thus the load vector is given by, 

{?} = {0 0g(a:,t/) 0 0 0 0} (2.37) 

q(x,y) is distributed load acting on the plane of the plate. For uniform loading 

q{x, y) = go (2.38) 


and for sinusoidally varying load. 


/ . . nx . Tty 

q[x,y) = gosm — sm — 
a 0 


Substituting for q(x,y) in expression for V in Eqn.2.36, we get 

. /NN \ 


(2.39) 


(2.40) 


where, 

{F}<‘> = y‘'’[iV]('l {«}<'> (2.41) 

F^®^ is the elemental load vector. It is obtained by assembling the nodal sub-vectors. 

The elemental load vector in terms of sub- vector can be written as, 

NN 

{n"’ = E{9i} P-«) 

t=l 

Where gi, is given by 

{ft} = 'tl:w,w,det[m {ft} (2 «) 

p=l?=l 

Hence the load acting on node of the element is obtained by 

9*3 — j2f:W,W,det{JWi9i3 

p=:i 


(2.44) 



2.8 Total potential 


For a conservative system, the total potential is stationary at the equilibrium configuration. 
If U is the strain energy or the internal potential developed and V is the potential of the 
system due external workdone, the total potential of the system is given by 


5{U + y) = 0 


5 n =0 


(2.45) 


Hence, 


JJ = C/ + y = Constant 


(2.46) 


Further we can write, 


5 [KY = 0 

q{ci}'iifiw+{<i}'{F}''>}=o 


( 2 . 47 ) 


where {d} is the global displacement vector, [K] is the global stiffness matrix and {F} is 
the global load vector. These are defined as 

NE 

Id} = E {<'}'•’ 


NE 

EM"’ 

e=l 


{F} = -Em' 


(2.48) 


Equation.2.47 may be rewritten as, 


S {d}' ((if) {d} + {F}) = 0 


(2.49) 
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The variational Eqn.2.49 leads to a set of linear equations of the form, 


[;i1W + {n = o 


(2.50) 


The above Equation 2.50 is the governing equation of flexure for laminated composite 
plates. This equation is solved incorporating appropriate boundary conditions to obtain the 
displacement vector. The formulation presented herein was implemented in HP9000 series 
machine. The governing equation has been solved using system(NAG) library routines. 



Chapter 3 


Solution approach - Monte Carlo 
technique 


The problen. of analysis of composite structures with variation in material properties usmg 
Monte Carlo approach can be broken down into the following steps t 


1. Random input sample generation. 

2. System analysis. 

3. Post processing of response behavior. 

These steps are represented by a block diagram in Fig.3.1 


3.1 Input sample generation 

3.1.1 Number generator 


A set of random numbers are generated of a given sample size, mean and standard de- 
viation. These numbers are then manipulated in sud. a way that these have a Gausstan 
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OUTPUT CHARACTERISTICS 


Figure 3.1: Lay-out of the code developed 
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3.1.2 Preprocessing of raw data 

To get an input sample with desired characteristics, the pseudo random numbers generated 
by G05FDF are preprocessed before feeding them as input to the FEM analyser. The 
samples are checked for shift in the input mean(/iinput) and standard deviation(a,„p„t) 
which are usually present in a small amount, depending on the arithmetic of the computer 
employed. 

1. Removal of shift in SD 

The standard deviation cr of a raw data can be obtained from the relation for the 
variance 

(S?_i one can write 

Each number is multiplied by the factor tsn to get a sealed sample with Odesired 
the standard deviation. 

2. Removal of shift in the mean 

Let Xi (z = 1, 2, . . . , n) be a set of pseudo random numbers, 
then, E"_iXi/n = ji 

where, /z is the mean of the generated numbers. 

After the sample is checked for the error in standard deviation, the sample with 
fJ-deaired is obtained by offsetting each of the Xj by a constant e, where emean=P- — 
/^desired- TFc Conditions for offsetting are ; 

if e < 0 the sample with Udesired is obtained by adding e to each number, in the case 
where e > 0 the numbers are reduced by e. 

The preprocessing is achieved by implementing the steps in succession. For this the 
statistical averaging subroutine (NAG-GOIAAF) is combined with the random number 
generator (NAG-G05FDF). These samples are subsequently used as the input to the finite 
element analyser to study the plate deflection behavior. 
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3.2 System analysis 

The governing equation of the bending of anisotropic laminated composites has been solved 
with the help of the finite element technique. A 7-DOF C° continuity model with 9-noded 
Lj^rangian element is employed here. The presented analyser have been organized into 
following modules. 

1. Mesh Generation Module. 

2. Constituent Matrix and Force Vector Module. 

3. The Solver. 

The plate is discretized into various elements. The mesh generation module generates 
the nodal co-ordinates for each element. Subsequent to this it generates the connectiv- 
ity matrix and identifies the elements and their respective nodes on the boundary. With 
this mesh generation scheme full plate and quarter plate analysis can be performed. The 
boundary conditions are appropriately applied for both the cases. The code also can gen- 
erate co-ordinates, connectivity and apply boundary conditions for plate with any shape 
of cut-out. The module generates a uniform mesh around the cut-out with quadrilateral 
elements. The mesh for the plate with cut-out is generated only for quarter plate. The rou- 
tine has the capacity to generate mesh for arbitrarily shaped cut-outs, provided the shape 
is represented as polar function f{r,d). As the thickness (h), individual ply-orientation 
and the material properties are given, a general constituent matrix is developed. With 
simple changes in the input data isotropic, orthotropic and anisotropic constituent matri- 
ces can be developed. This along with the element co-ordinates go to form the element 
stiffness matrices. Element stiffnesses are assembled with the help of connectivity matrix 
and boundary conditions are appropriately applied to the boundary elements. With the 
information about the nature of loading the element force vector is generated and assem- 
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bled. The assembled global stiffness matrix and the global force vector constitute a set of 
linear equations. NAG F04ATF solver is used to get the respective deflections depending 
on the type of the plate analyzed. The deflections obtained from the developed FE model 
has been validated against results of standard HSDT models [13, 14]. 

3.3 Post processing of the response behavior 

The sample of the deflection response of the plate is obtained from the system FEM anal- 
yser. Its statistical information is extracted by feeding the obtained data sample to the 
routine GOIAAF. This routine gives the mean and the standard deviation of the fed data. 

The above scheme is used for obtaining the displacement characteristics with the ba- 
sic material properties varying individually as well as simultaneously for plates with and 
without cut-outs. These are presented graphically in Chapter 4. 
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Chapter 4 


Results and Discussion 

The approach presented in the previous chapter is used to obtain the deflection statistics 
of rectangular plates. This chapter is organized in two sections. In the first section the 
behavior of rectangular plate with different support conditions has been studied for change 
in material property variation over a range and in the next section the study is conducted for 
r«:tangular plates with cut-outs. Initially the efltect of randomness in individual material 
property has been studied to identify the most critical parameter. In the latter section a 
problem of practical importance is dealt with assuming all the material properties to vary 
randomly and simultaneously. 

• Samples are generated with the standard deviation varying from 2 percent to a of 20 
percent of their mean values with sample-size 1000 and 10 such samples are used for 
each study. 

• The material properties considered to be random are Ei.EuViuGu and Gtz- 

Results have been obtained for T300/52208 graphite/epoxy composite [6] with mean values 
as given below ; 
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E, = 145.0 GPa Gu = 4.5 GPa 

Et = 10.71 GPa Gu = 4.5 GPa 

I'll = 0.31 Gu = 3.5906 GPa 

4.1 Deflection characteristics of rectangular plates 

In this section we study the effect of dispersion in individual input random variables while 
the characteristics of the other input variables are kept constant on plates with different 
aspect-ratios (AR), boundary conditions, ply-orientations and thickness-ratios (a/h). Here 
the variations are considered on a SS square antisymmetric laminate with ply-orientation 
[45°/ — 45°/45°/ — 45°] and a/h of 10. In a given study except for the parameter selected all 
the other parameters are kept constant. The ARs chosen for the study are 2, 3 and 4. Two 
classical boundary conditions have been considered namely all edges SS and all edges CC. 
Results are obtained for 4-layered antisymmetric laminates with ply-orientation 15°, 30° 
and 45° and 0°. Three a/hs of 4(thick), 8 and lOO(thin), are considered for the plate. The 
curves for ajh = 100 are multiplied by 10“^ in the case of Eu Et and uu characteristics, 
and for the characteristics of Gu and Gu the factor is 10~\ 

Plate is uniformly loaded with q = lOOOiVm"^ and displacement refers to central dis- 
placement of the plate. 

The Fig. 4.1 through Fig. 4.5 shows the response of the plate to each input RV. Each 
figure represents four different studies A, B, C and D. where A, B, C and D presents the 
effect of the SD of input RV on the SD of the central deflection of the plates with different 
boundary conditions, ARs, a//is and ply -orientations. 



Effect of Longitudinal Elastic Modulus Ei 


Fig.4.1 gives the characteristics of the plate response with different standard deviation of Ei. 
The response shows a linear behavior for small cr^,. However, with increase in deviation 
the response becomes somewhat non-linear. The non-linearity is stronger for some ply- 
orientations compared to the effects of support conditions, ARs and a/h. The dispersion 
in the displacement characteristics of the plate is found to be maximum for a/h = 100 
(thin plate). The plate shows least sensitiveness to input dispersion for the AR = 4. There 
is a marked increase in sensitiveness of the plate response for AR=2 compared to the 
AR=3 and 4. The plate with all edges CC is more affected by change in input randomness 
compared to a plate with all edges SS. From the study of the effect of a/h it is found that 
variation in the characteristics is most prominent for thin, and least for thick plates. The 
0° symmetric laminate is affected more and a 45° anti-symmetric laminate gets affected 
the least due to change in However, the overall difference between the two is small 
indicating low sensitivity to ply-orientation. 
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Figure 4.1; Plate deflection characteristic with as, 
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Effect of transverse elastic modulus Et 

The characteristics of the system with Et as input random variable is presented in Fig.4.2. 
The general behavior is almost linear except for the plates with different ply-orientation, 
but unlike the case of Ei, here its slope shows a decreasing trend. The plate with AR=4 
is least affected by change in the input randomness while the plate with ajh— 100 shows 
maximum effect. As observed in the case of Ei there is again a marked higher sensitiveness 
for plates with AR=2 as compared to plates with AR=3 and 4. All edges simply-supported 
plate is more affected by dispersion in Et than all edges CC plate. The effect of the a/ his 
similar to that observed for Ei. The plate with 15*’ orientation is affected maximum and 
indicates higher sensitivity compared to the 0** symmetric case. The slope of the response 
characteristics are lower for the other two angles, reducing from 30** to 45*’ anti-symmetric 
case. 

Effect of major Poisson’s ratio uu 

The sensitivity of mid point displacement characteristics of the plate for variation in 
is given in Fig.4.3. The general nature of the curves are mostly linear. The least affected 
are the plates with a/h = 4 and plate with AR=4. On the other side the plates with 
ply-orientation 0** and ajh = 100 get affected the most. The effect of change in the AR 
shows a similar trend as discussed before for Ei and Et- There is a reversal of trend 
for the two support conditions compared to the effect of Ei and Et as the displacement 
characteristics of SS case is affected more than that of the CC edge one. As the ply 
orientation angle increases the response slope reduces, the reduction being stronger for 
small angles, indicating decreasing trend with increase in the ply angle. 
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Figure 4.2: Plate deflection characteristic with 
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Effect of rigidity modulii Gu and Gtz 

Fig. 4.4 shows the behavior the system for input random variables Gu and Gtz- The 
characteristics are mostly linear. The 0° symmetric laminate and plate with a/h = 100 is 
affected maximum and the plate with AR=4 is affected the least for change in dispersion 
in Git and Gtz- The trend of all the characteristics, in general are similar to that of the 
Pit case . The response for a/h of 4 and 8 are close and much less compared to the case of 

a/h = 100. 

Effect of rigidity modulus Giz 

The system behavior for input random variable G„ is given by Fig. 4.5. Here too the 
characteristics are mostly linear and plate with a/h = 100 is affected the most and plate 
with AR of i is least affected. The response curves for both the boundary conditions almost 
coincide with each other, showing that behavior is not affected with the two boundary 
conditions considered. The behavior of all the other characteristics, in general resemble 

the Git case. 
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4.1.1 Response sensitivity 

In this section, a comparative study between different plate configurations is carried out to 
highlight the quantitative sensitivity of the system to randomness in individual material 
properties. The transverse loading is assumed as uniformly distributed 1000iVm~^ and 
unless specified the ply-orientation is [45°/ - 45°/45°/ - 45°]. Four plates are used for the 
sensitivity study : 

1. A square plate, with ajh of 10, CC at all edges. 

2. A rectangular plate with AR=2 and a/h=10, with SS edges. 

3. Square and moderately thick plate a/h=8, with SS edge conditions. 

4. A square plate, with SS edges, a/h = 10 and ply-orientation [15°/ - 15°/15°/ — 15°] 
and a/h = 10.. 

Fig.4.6 shows the deflection characteristics for the four plates for variations in the standard 
deviation of the material properties Ei, Et, uiu Gu and Gu . The material properties are 
taken to vary only one at a time, other being held constant at their mean values. It is seen 
that for the cases studied here the plate is most sensitive to the randomness in the longi- 
tudinal elastic modulus(jBi), the next most sensitive parameter being the rigidity modulus 
Git sod Gtz, The remaining material parameters do affect the displacement characteristics 
but the effect is less than half the effect of the longitudinal modulus. 
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Figure 4.6: Response sensitivity of the comp<»ite plate to material properties 
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4.1.2 Normalised characteristics of composite plate 

Effect of variation in standard deviation of the individual material properties, keeping 
the rest all constant, has been presented in detail in the preceding sections. Prom the 
applications point of view it is appropriate to consider the case where ail the properties 
vary simultaneously. Fig.4.7 presents the normalised characteristics of the plate, with all 
the material properties varying simultaneously each assuming same value for the ratio of 
its standard deviation to mean. 

The displacement standard deviation is normalised with the mean displacement of the 
plate obtained with all the material properties having their mean values. Fig. 4.7 shows the 
normalised behavior of the plate. It is observed that the proposed normalisation scheme 
when used, the different characteristic curves for different ARs happen to fall on each other, 
this is also true for the characteristic curves for ajhs and boundary conditions. The plate 
with [45°/ — 45°/45°/ — 45°] ply-orientation is least affected by the input randomness, 
and the 0° symmetric laminate has higher response slope. However, the dispersion in 
response characteristics for different ply-orientations is small. Fig. 4.8 is constructed by 
superimposing all the normalised curves. One can observe that the curves form a band 
which is non-linear in nature. The non-linearity enter into the characteristics beyond 10 
percent of the input a. The maximum dispersion in the characteristics is 1 percent and 
the curves lie in between a maximum ffoutput of 23 to 24 percent for the normalised Omput 
of 20 percent. 

The interesting finding has raised queries whether the normalised characteristics of all 
plates behave in a similar way. Hence, further studies are pursued in this direction for 
plates with cut-outs and a different loading pattern. 
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4.7: Normalised reponse characteristics of rectangul^ plates 
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Figure 4.8: Superimposed plate response characteristics 
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4.2 Deflection of plates with cut-outs 

In Aerospace structure applications weight reduction is an important applications design 
criteria. Besides being used due to functional requirements, cut-outs also have relevance 
to the problem of weight reduction. More and more structures are used with cut-outs. 
This motivated the study of the behavior of composite plates with cut-outs. The analysis 
approach developed can handle different shapes of cut-outs. Here, for the present study, 
circular cut-outs have been used in rectangular plates with plate length to cut-out diameter 
ratios or cut-out ratios {l/d) having values 2,4, 6, 8 and 10. A standard plate is taken 
with AR = 1, a/h = 10, ply-orientation [45®/ — 45®/45®/ — 45®]. A sinusoidally varying 
distributed load of amplitude lOOOA’ is used for this analysis. The displacement of the 
point located at the intersection of rim of the cut-out and the major axis of symmetry is 
monitored and taken for the study. 

The normalisation scheme of the previous section is used here for the study. Plots are 
developed to bring out the effects of boundary conditions, ARs, a/h and ply-orientation 
for the five l/d ratios mentioned. Figs. 4.9. to 4.13 present normalised displacement char- 
acteristics of the plates with different l/d ratios. The use of this scheme brings forward 
anticipated interesting results. In these figures, it is observed that irrespective of the dif- 
ferent configurations of the plate, such as different AR, boundary conditions, ajh and 
ply-orientations, their normalised displacement characteristic curves are very close, almost 
coincide with each other. 
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Figure 4.10: Characteristics of the composite plate with cut-out, Ijd = 4 


a* 







44 







45 



Figure 4.12: Characteristics of the composite plate with cut-out, l/d = E 
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This aspect is further illustrated by the Fig. 4.14 , obtained by superposing all the Fig.4.9 
through Fig.4.13. A marker line y = x is included to show the deviation of these curves 
from linearity. The characteristic curves group themselves into two bands. The thicker 
band representing the behavior of these plates with SS edge condition and the band with 
maximum deviation is the characteristics of the plates with cut-out l/d ratios ranging from 
2 to 10 and all edges CC. It is observed that normalised response characteristics is affected 
by different boundary conditions and the nature of the curves are linear up to 10 percent 
of the variation in the input characteristics, beyond this the the curves show a distinct 
non-linearity the deviations being, 23 percent for the SS plate and 24 percent for the CC 
plate for the input SD of 20 percent. 

It can be inferred from these figures that the effect of variation in the randomness in 
the material property on different plates with different AR, a/h, support condition and 
fiber orientation may be presented in the normalised form in one similar way. In other 
words, the curve may be thought as a design curve which can predict the normalised char- 
acteristics for all the configurations mentioned. It is understood that from the normalised 
characteristics we can always obtain the absolute characteristics of any plate configuration. 
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Chapter 5 


Conclusions and scope for future 
work 


In this thesis, a statistical study of the bending behavior of composite laminated rectangular 
plates with different configurations has been carried out using finite element method. From 
this study the following main conclusions can be drawn. 

1. The longitudinal modulus Ei and the rigidity modulus Gh are the most critical of 
the material properties. They considerably influence the deflection characteristics of 
the composite plate. 

2. The thin plate (a/h = 100) is most sensitive and the plate with AR=4 is least 
sensitive to the variation in input material property characteristics. 

3. The behavior of the plate is found to be linear for low dispersion in material proper- 
ties. For larger dispersion, with {^) becoming greater than 10 percent, the behavior 
turns non-linear. 

4. The normalised characteristics of the plates is affected by the end conditions. The 
ply-orientation also affects the characteristics but its eflfect is quite small compared 
to the effect of end conditions. 
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5. With 20 percent of input SD the change in SD of normalised deflection characteristics 
is 23 percent for the SS plate and 24 percent for the CC plate. 

6. The normalised characteristic curves of the plates with different ARs, boundary con- 
ditions, a/h and ply-orientations behave in a similar way and all of the curves fall 
within a tolerance of 1.5 percent of Ooutput- 

7. With one normalised behavioral curve, the displacement characteristics for a family 
of plates can be predicted. The actual displacement characteristics for these can be 
derived as the normalised factor is known. 


Scope for future work 

Having worked with this interesting problem for long the author would like to suggest the 
following for future work: 

1. The effect of dispersion in the random material properties on the stresses and dy- 
namics of the plate can be studied. 

2. The random material properties were taken to have Normal distribution over lots 
and was assumed not to vary over the domain of the plate. The problem may be 
pursued with the material properties having different probability distributions over 
the domain of the plate as well as varying in lots. 

3. The above problem can be pursued with random loading, plate geometry being circu- 
lar, elliptical etc. The boundary condition may be fixed, SS, free and in combination 
at different edges. The parameters a/h, fiber orientation, etc. may be taken as 
random in nature. 



References 


[1] S. Salim, “Analysis of composite plates with random material properties.,” Phd-Thesis 
Dept of Aerospace Engg., IIT Kanpur., April 1995. 

[2] R. A. Ibrahim, “Structural dynamics with parameter uncertainities,” Trans., ASME; 
Applied mechanics review, vol. 40, no. 3, pp. 309-328, 1987. 

[3] S. Nakagiri, H. Takabatake, and S. Tani, “Uncertain eigen value analysis of composite 
laminated plates by the sfem,” Composite Structures, vol. 109, pp. 9-12, 1987. 

[4] A. W. Leissa and A. F. Martin, “Vibration and buckling of rectangular composite 
plates with variable fiber spacing,” Composite Structures, vol. 14, pp. 339-357, 1990. 

[5] S. P. Englestad and J. N. Reddy, “Probabilistic methods for the analysis of metal- 
matrix composite,” Composite Science and Technology, vol. 50, pp. 91-107, 1994. 

[6] G. V. Vinckenroy and W. P. de Wilde, “The use of monte carlo techniques in sfem 
for determination of the structural behavior of composites,” Composite Structures, 
vol. 32, pp. 247-253, 1995. 

[7] S. Salim, D. Yadav, and N. G. R. Iyengar, “Analysis of composite plates with random 
material characteristics.,” Mechanics Research Communication, vol. 20, no. 5, pp. 405- 
414, 1993. 

dNTRAL LiBRAR' 

I !. T.. KANPIM* 


51 


[8] S. Salim, D. Yadav, and N. G. R. Iyengar, “Deflection of composite plates with ran- 
dom material characteristics.,” Proceedings of Int. Sym. on Aerospace Science and 
Engineering., pp. 236-239, Dec 12-15 1992. 

[9] S. Salim, D. Yadav, and N. G. R. Iyengar, “Free vibration of composite plates with 
randomness in material properties.,” Proceedings of the Fifth Inti. Conf. on Recent 
advances in Str. Dynamics, pp. 814-823, July 18-21 1994. 

[10] R C. Yang, C. H. Norris, and Y. Stavsky, “Elastic wave propagation in heterogenous 
plates,” Int. J. Sol. Struct, vol. 2, pp. 665-684, 1966. 

[11] K. H. Lo, R. M. Christenson, and E. M. Wu, “A high-order theory of plate deforma- 
tion,” Trans. ASME, J. Appl. Mech., vol. 44, pp. 663-668, 1977. 

[12] J. N. Reddy, “A simple higher order theory for laminated composite plates,” Trans. 
ASME, J. Appl. Mech., vol. 51, pp. 745-752, 1984. 

[13] C. A. Shankara, “Analysis of composite laminated plates subjected to mechanical, 
hygral and thermal loadings.,” Phd-Thesis Dept, of Aerospace Engg., IIT Kanpur., 
pp. 43-81, September 1993. 

[14] J. N. Reddy, “A refined non linear theory of plates with transverse shear deformation,” 
Int. J. Sol. Struct., vol. 20, pp. 881-896, 1984. 



